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CREEP BENDING AND BUCKLING OF LINEARLY 
VISCOELASTIC COLUMNS 
By Joseph. Kempner 

SUMMARY 


The general dynamic equation of creep bending of a beam loaded lat- 
erally and axially was derived for a linearly viscoelastic material whose 
mechanical properties can be characterized by four parameters . The mate- 
rial can exhibit instantaneous and retarded elasticity as well as pure 
flow . 


The equation derived was -used to obtain the creep bending deflection 
of a beam in pure bending and of a column with initial sinusoidal devia- 
tion from straightness. As expected, the ratio of the creep deflections 
of the beam in pure bending and the deflections of a corresponding purely 
elastic structure is identical to the ratio of the creep strain and the 
corresponding elastic strain of a bar under simple tension or compression. 

The results of the analysis of the creep deflection of the column 
showed that the deflections increase continuously with time anrl become 
infinitely large only when the loading time is correspondingly large . 
However, large deflections are obtained in reasonably short periods of 
time if the applied load Is near to the Euler load of the column. Th<* 
deflection-time curves obtained from a numerical example are of the same 
type as those determined by experiment with aluminum col umn s . 


INTRODUCTION 


The problem of the determination of the behavior of beams pmrl columns 
under conditions conducive to creep deformation has been given attention 
only in the past several years (see, for example, refs. 1 to 8). This 
problem is becoming increasingly important because of the high tempera- 
tures at which high-speed missiles and aircraft as well as modern power 
plants operate. It is the purpose of the present report to investigate 
analytically the nature of the creep deformation of beams formed from 
materials which can be assumed to exhibit linearly viscoelastic properties 
either at room or elevated temperature s . 
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The fundamental behavior of ideal viscoelastic materials when, sub- 
jected to various types of loading has been discussed rather fully in 
the literature (see, for example, refs. 2 and 9)* In this work the vis- 
coelastic material is represented by a system of linearly elastic ele- 
ments (such as Hookean springs) and linearly viscous elements (Newtonian 
dashpots) connected in various series and parallel configurations (see 
figs. l(a) and 1(b)) . An important class of problems which may be ana- 
lyzed with the aid of models of the aforementioned type consists of those 
involving the effect of viscoelastic creep on the response of beams and 
columns to various types of loads . Freudenthal (refs . 1 and 2) has dis- 
cussed such problems in which the structural material is considered to 
be represented by a single Maxwell element consisting of a linear spring 
connected in series with a linear dashpot; hence the material properties 
are defined by two parameters (fig. 1(c)) . However, it is shown in ref- 
erence 9 that in order to approximate real materials, such as amorphous 
linear polymers, by a model consisting of springs and dashpots several 
such elements must be incorporated in the model. Thus it is the purpose 
of the present report to extend the application of the concept of repre- 
senting the mechanical properties of viscoelastic beams and columns by 
networks of springs and dashpots to models characterized by more than 
two parameters . In particular, the general dynamic equation of motion 
of a bent beam column is derived for a material represented by a Maxwell 
element connected in series with a Kelvin or Voigt element where the 
latter element consists of a spring and dashpot in parallel (fig. l(d)), 
and hence the complete model involves four parameters (fig. 2). Such a 
model Is capable of exhibiting instantaneous and retarded elastic response 
as well as pure viscous flow. 

The author Is indebted to Professor N. J. Hoff for his guidance and 
criticism, to Mr. S. A, Patel for his assistance in carrying out calcula- 
tions, and to the National Advisory Committee for Aeronautics for sponsor- 
ship of the research reported in this paper . 


SYMBOLS 


A 

A,B,Ci,C2 

A',B' 


cross-sectional area 
constants of integration 

constants 


a 


Ei,E 2 


amplitude of function representing initial deviation from 
straightness 

moduli of elasticity 


G 


function of time 
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I 

k = 1 + 

L 

M 

M 

Mo 

m 

P 

P E 

P 

T 

t 

w 

w elastic 


moment of inertia of cross-sectional area of beam or 
column with respect to a principal centro idal axis 

( T a/ T i) + ( e i/ e 2) 

length of beam or column 
bending moment 

bending moment due to lateral and inertial loads 
constant bending moment 
mass per unit length 
constants 

axially compressive end load 
Euler load, 

applied lateral load per unit length 

function of time 

time 

deflection due to loads 
elastic deflection due to loads 


w i 


function representing initial deviation from straightness 


x 

y 

z 

a 

€ 

e elastic 


coordinate along centro idal axis of beam or column 

principal axis of cross-sectional area of beam or column 

distance from xy-plane 

ratio of end load and Euler load 

strain, positive in tension 

elastic b train 


strain at z = 0 
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e l' € 2> e la> € lb 


strain contributions 


\ l *^2 

P 

a 

ff 2a * a 2b 


coefficients of viscosity- 

radius of curvature 
stress, positive in tension 
stresses 


t^,T 2 relaxation and retardation times, respectively 

(■)•&()/* t 

( ) x H &( )/dx 


DERIVATION OF DIFFERENTIAL EQUATION OF BENDING 
Basic Theory 


The mechanical behavior of viscoelastic materials, such as amorphous 
linear polymers, under conditions of small stresses and strain, has been 
discussed in great detail by Alfrey in reference 9 • It iB sho-wn in thiB 
reference that an analogy exists between the mechanical behavior of such 
materials and the mechanical behavior of Hookean springs (linear elasticity) 
and Newtonian dashpots (linear viscosity) connected in various series and 
parallel networks (fig. l) . If two or more elements are coupled in series 
producing a configuration such as the Maxwell model of figure l(c) each 
of the elements must support the total load applied to the model, and 
the total extension of the system is the sum of the individual extensions 
of the elements involved. On the other hand, if elements are coupled in 
parallel (such being the case with the Kelvin or Voigt model of fig. 1(d)), 
the total load supported by the model is equal to the sum of the loads 
carried by the individual elements, whereas the total extension of the 
model is equal to the extension of each of its components . The Maxwell 
model is capable of exhibiting instantaneous elastic response and viscous 
flow when subjected to arbitrary loading conditions . The former action 
produces recoverable deformation, whereas the latter action produces per- 
manent deformation. The phenomenon of stress relaxation can be repre- 
sented with the aid of such a Maxwell model. The Voigt model represents 
a retarded or time -dependent elastic response to applied loads and is 
capable of exhibiting elastic aftereffects . It is shown in reference 9 
that the mechanical behavior of combinations of Voigt models in series, 
such as that shown in figure l(a) where the first and last elements may 
be considered as degenerate Voigt elements, or of combinations of Maxwell 
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elements in parallel (fig. 1 (b)) is analogous to the mechanical behavior 
of amorphous linear polymers under conditions of small stresses and 
strains . In order to determine the pertinent elasticity and viscosity 
coefficients to be used to complete the choice of a model representing 
a given material, suitable creep or relaxation tests must be undert ake n - 
The various coefficients are then determined empirically. Once these 
coefficients are known as f mictions of temperature (and they will be func- 
tions only of temperature if the material being represented is visco- 
elastic and the model chosen is of sufficiently complete structure) , then 
the behavior under arbitrary conditions of loading and temperature of 
any isotropic, homogeneous structural member formed from the viscoe las tic 
material can be described analytically. Although only those structural 
members which can be considered to be uniaxially stressed are considered 
herein, a generalized three -dimensional theory exists and is described 
in reference 9 . However, the theory presented therein is valid only for 
models characterized by two parameters . 

In the ensuing analysis the behavior of viscoelastic beams nnrl col- 
umns is studied by assuming that the material of construction can be 
represented adequately by a network consisting of a spring in series with 
a Voigt element and a dashpot (see fig. 2) . The model chosen can be con- 
sidered as representing a first approximation to a more general model 
(see fig. 1 (a)) and is capable of exhibiting instantaneous and retarded 
elastic response as well as viscous flow. Any model of simpler nature 
than that chosen is in general incapable of responding to loads in a 
manner analogous to known viscoelastic materials . In the present study 
the structures considered are assumed to be in a state of uniform con- 
stant temperature and hence the elasticity and viscosity coefficients are 
constants . 


STRESS -STRAIN-TIME LAW 


* 


With the use of the principles governing the behavior of springs and 
dashpots coupled in series or parallel stated earlier, the general stress- 
strain-time relationship governing the behavior of the four -parameter model 
of figure 2 can be determined. If at any time t the stress a is 
applied to the model as shown in figure 2 , where a is positive in ten- 
sion, the strain contribution e-^ of the elastic element (la), having 

an elastic modulus E]_, is given by Hooke's law as 


e la = */ E l (1) 

The relation between the applied stress and corresponding strain in the 
viscous element (lb) is given by Newton's law of viscosity as 


4 lb = ff / A l 


( 2 ) 
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in which is a coefficient of viscosity and the dot above the strain 

indicates differentiation with respect to time. Thus the combined strain 
contribution of elements (la) and (lb) is related to the applied 

stress as follows : 

El = + ( ff Al) ( 5 ) 

Equation ( 3 ) is the well-known stress -strain-time relationship for a 
Maxwell model. 

The strain contribution £2 of elements (2a) and (2b), respectively, 
is related to the applied stresses a 2a a 2b> w ^ iere "kk® sum these 

stresses is the applied stress c, by ■ expressions similar to equations (l) 
and (2) . Thus 


e 2 = a 2a/ E 2 

(4) 

e 2 = °2b/ A 2 

(5) 


in which E 2 is the modiolus of elasticity and Ag the coefficient of 
viscosity of the elastic and viscous elements (2a) and (2b), respectively. 

Hence the strain e 2 is related to the applied stress by the 
equation 

A 2 € 2 E 2 e 2 = ^ ( ^) 

Equation (6) is the defining equation of a Kelvin or Voigt model. 

The total strain in the model of figure 2 is 

e = e 1 + € 2 (7) 

in which and e 2 are defined by equations ( 3 ) and (6), respectively. 

Elimination of £-j_ and e 2 from equations ( 3 ), ( 6 ), and ( 7 ) leads to 
the followin g stress-strain-time relation governing the mechanical behavior 
of the four-parameter model considered: 

a + (k/r 2 )a + (i/t^o = E i ( 8 ) 

in which t^_ = A^yE-^, the relaxation time of a Maxwell element with elas- 
ticity and viscosity coefficients Ej_ and A]_> respectively; t 2 = ^ 2 / E 2 > 
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the retardation time of a Voigt element with elasticity and viscosity 
coefficients E 2 and Ag, respectively; and k = 1 + + ^E-j^E 2 j. 

Equation (8) is thus the uniaxial stress-strain -time relationship 
for a viscoelastic material whose response to loading at a constant tem- 
perature is analogous to the behavior of the model shown in figure 2. 

The determination of the strain-time relation for a member under uniform 
tensile or compressive stress is facilitated by the direct application 
of the component equations ( 3 ) and (6) rather than the general equation (8). 
Thus equations ( 3 ) and (6) together with equation ( 7 ) and the condition 
that at t = 0 € is completely elastic and equal to result in the 

following strain-time relation (see ref. 9) : 


e / e elastic = ( E l/ E 2) (p- " e ^ 2 ) + ( t / T l) + 1 

in which the elastic strain 6 e lastic = a j^ 1* ^kis relation is plotted 

in figure 3 for values of the elasticity and viscosity coefficients used 
in subsequent calculations . 


Strain-Deflection Relation 

Consider a beam loaded in a manner such that the bending moment 
vector M acting on any cross section coincides with a principal cen- 
troidal axis such as the y-axis in figure b. With the assumption that 
cross sections of the undeformed beam remain plane and normal to the cen- 
troidal axis after bending, each element of the beam deforms under the 
influence of the bending moment and an axial load as indicated in fig- 
ure 5. In this figure x is the centro idal axis of the beam; z is 
the distance from the xy-plane to any other parallel plane; e Q and e 

are the strains, shown positive in tension, at the xy-plane and at a 
distance z from this plane, respectively; and p is the radius of 
curvature of the centro idal axis . Thus the strain corresponding to any 
z is 


6 = (1 + € 0 )(z/p) + e Q ( 9 ) 

If the strains are considered small, equation ( 9 ) becomes 

e = (z/p) + € q (10) 

Since the strain e is a linear function of z, the stress is also a 
linear function of z (see eq. (8)). Hence, the stress distribution 
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across a section can be decomposed into a uniform stress and a l ine arly * 

distributed stress which vanishes at the xy-plane . If a constant axially 
compressive load P is applied to the beam of cross-sectional area A, 
the uniform stress (and hence the stress at z = 0) is* -P/A. Therefore, 
e 0 , which is caused by the constant stress -P/A, can be determined from 

equation (8a), in which ^elastic ■'■ s now 'taken a® -P/ AEq . 

For small displacements w of the centroidal axis 

l/p = -Wjrjr (H) 

where a subscript x denotes one differentiation with respect to x. 

Thus for small strain and displacements 

e = - zw xx + e o C 12 ) 


Differential Equation of Bending of a Beam Column 

With the aid of equations (8) and (10), the differential equation 
of bending is readily derived. Multiplication of equation (8) by z dA 
and integration of the resulting equation over the area A yield 


J ffz dA + (k/r 2 ) J ^A + (l/ T l T 2) J 


crz cLA = 


Eq J ez cLA + ^ 1 /^ 2 )^ ^ 


ez dA 


( 15 ) 


The moment M is related to the stress by the relation 


M = / cz dA 

'A 




Substitution of equations (10) and (l4) into equation ( 13 ) together with 


z £ dA = I 



and 
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where I is the moment of inertia of the cross-sectional area 
respect to the centroidal principal axis y (fig. ^4-) , yields 


A with 


E X I 


d 2 (l/p)y4t 2 + (l/T 2 )d(l/p)y&tJ = M + (k/r 2 )M + (V T i T 2) M (l?) 


Equation (15) is the general equation of bending of a beam of viscoelastic 
material whose mechanical analog is the model shown in figure 2. For 
small displacements the curvature l/p is defined by equation (ll) and 
hence equation (15) reduces to 


-Eilj^xx + (i/t 2 )^ “ ** + (V T 2)M + (i/ti t 2 )m 

The bending moment can be expressed as . 

M = P(w + w^) + M(x,t) 


(16) 


(17) 


in which w^ represents an initial deviation from straightness of the 
original beam axis and M is the moment contribution due to applied 
lateral and inertial loads. Substitution of M from equation (17) into 
equation (l6) yields 


El 1 j^xx + (V T 2 )^ + + (V t 2 )w + (i/tiT 2 )^ = 

-( P / T l T 2 ) w i “ “ " ( k / T 2 )^ ' (V T 1 T 2 )“ 


(18) 


Furthermore, if m is the mass per -unit length of the beam and P is 
an applied load per unit length of the beam, where m may vary with x 
and p with x and t, then 

Mxx - -p + mw (19) 


Hence equation (18) becomes 



(k/r 2 )p + (l/T^ivJp (20) 

Equation (20) is the general dynamic equation of bending with small deflec- 
tion of a four -parameter viscoelastic material. If t 2 is considered 

infinitely large, that is, the Voigt element of figure 2 is frozen in, 
then equation (20) reduces to the dynamic equation of bending of a beam 
of Maxwell material (see ref. 10) . 
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INITIAL CONDITIONS 


Regardless of the exact nature of any particular beam problem, cer- 
tain in itial conditions which must be imposed on all solutions of equa- 
tion (20) are assumed in the present application of the viscoelastic 
model of figure 2. These conditions are that all external loads are 
applied instantaneously and hence the initial response is purely elastic 
(see ref. 9) • Thus at t = 0 

6 = e la = °/ E l 

and 

e lb = e 2 “ 0 


( 21 ) 


( 22 ) 


Since € * e^ a + + e 2> equations ( 21 ) and (22) yield only two independ- 

ent relations. Differentiation of equation ( 7 ) with respect to t and 
subsequent substitution of equations ( 3 ) and (6) into the resulting equa- 
tion yield 

e 2 = ( T 2/ E l)<* + |(> ‘ ^l] 0 ~ t 2^ (25) 


Introduction of equation (12) 
relation for ^ 1 

e 2 = ( T 2/ E l) 0 + 


into equation ( 23 ) results in the following 

[jk - lyE^jb + t 2 (ot xx - e D ) (24) 


Multiplication of equation (24) by z dA, integration over the area A, 
and introduction of equation (l4) into the relation obtained yield for 
€£=0 at t = 0 



M + t 2 Iw xx = 0 


(23) 


Hence equations (21) and ( 25 ) are initial conditions which must be imposed 
on all solutions of equation ( 16 ) or ( 20 ) . 


APPLICATIONS OF DIFFERENTIAL EQUATION OF BENDING 
Beam Under Pure Bending 


Equation ( 16 ) can be readily applied to the problem of creep due to 
pure bending of a viscoelastic beam with inertial effects neglected. If 
the applied constant bending moment is designated as Mq, then equation (l 6 ) 
becomes 
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Ellj^xx + (l/ T 2)^xx] = (l/ T l T 2) M o 


( 26 ) 


Integration, of this equation twice with respect to x and introduction 
of the conditions that, at x = 0 and L, y = w = 0 yield 

Eiljy + (l/T^ = (Mo / /2t 1 t 2 )(Lx - x 2 ) (27) 

in which L is the length of the beam. 

If w is written as 

w = G(t)(M 0 /2T 1 T 2 )(l/E 1 l)(Lx - x 2 ) (28) 

where G(t) is a function of time only, then equation ( 27 ) reduces to 

G+(i/t 2 )g=1 (29) 

The general solution of equation ( 29 ) is 

-(V T 2 ) 


G = C]e + r 2 (t + C 2 - t 2 ) 


(30) 


where and C 2 are constants of integration. Thus equation ( 28 ) 

becomes 


w = (M 0 /2T 1 T 2 )(l/E 1 l)(Lx - X 2 ) 


C^e 


-(V T a) 


+ T c (t + Co - T, 


0 


( 31 ) 


At t = 0, the initial deflections are purely elastic and hence 

" = ^elastic = (M 0 /2E 1 l)(Lx - x 2 ) (32) 

Also at t = 0, e 2 = 0 and hence equation ( 25 ) together with equation ( 31 ) 
yields 


C 1 = _ ( E l/ E 2) T l T 2 


(33) 


Application of equation (32) and substitution of equation (33) into equa- 
tion ( 31 ) result in the following value for C 2 : 
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Hence from equations ( 31) > ( 33) * and ( 3*0 


(mo/sb-li)^ - x 2 ) (35) 

In terms of the elastic displacement, equation (35) becomes 

w / w elastic = (Ei/E 2 )( 1 - e ^ 2 ) + { t / T l) + 1 (? 6 ) 

Equation (35) or ( 3 6 ) describes the pure-bending creep behavior of 
a viscoelastic beam whose material has a mechanical behavior under stress 
analogous to the behavior of the model of figure 2. Equation ( 36 ) is 
plotted in figure 3 end yields a creep curve identical to the creep curve 
for the same material under simple tensile or compressive stresses (see 
eq. ( 8 a)) . 



Column With Initial Curvature 

The behavior of a viscoelastic column whose centroidal axis initially 
deviates from a straight line can also be analyzed with the aid of equa- 
tion (20) . With inertial effects neglected equation (20) when applied 
to the simply supported uniform column of length L indicated in fig- 
ure 6 reduces to 



If the initial deviation from straightness is taken as 

= a sin(jtx/L) ( 38 ) 

where a is a constant, then the solution of equation (37) for the 
additional deflection w can be assumed in the form 

w = aT(t) sin(«x/L) (39) 

where T(t) is a function of time only. Upon substitution of equations ( 38 ) 
and ( 39 ) into equation ( 37 ) the following equation results for the func- 
tion T: 
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(V T i T a )(*/%)/£ - (,,o) 

in which the Euler load P E = n^Ej_l/L^. The general solution of equa- 
tion (40) is 


nit ngt 
T = Ae + Be - 1 


(hi) 


1/2 


in which A and B are constants of integration and 

n l - -[(V^a)^ - to )/( 1 - a D + I 1 / 21 !)! 1 - to)/(l - a]] 
{(Vvs^/l 1 - “J] 

n 2 - - te)/(l - 4] - f [[l/2T 2 )(l - ta)/(l - a]] 1 

j7 i / T i T 2) a /( i - <*- y| J 


)l/2 


(h2) 


where a = jEy'Pg : Since a must be less than unity in order for the 
column to be initially stable, n^ is real and positive and n 2 is 
real and negative . 

Upon substitution of equation (4l) into equation (39 ) > the deflec- 
tion w becomes 


/ n-jt n 2 t 
w = ( Ae + Be 

The initial conditions which must be 
eqs. (21) and ( 25 )) 

t = 0, w = Veias-fcic = 

t = 0, (Tjj/E^M+^k 


- lja sin(jtx/L) 
/ 

(43) 

imposed on equation (43) 

are (see 

ja,/ (l - <x)j a sin(itx/L) 

(44) 

“ i)/sJm + Tglw^ = 0 

( 45 ) 
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Since M = P (w + w-j^, equation (45) together with equation (43) yields 
at t = 0 

j*l - [( k " 1 )/ T 2] a /( 1 “ a )j A + |^ n 2 " [(* " ^./(l - a) j B = 0 (46) 

Similarly, equations (43) and (44) yield 

A + B = 1/(1 - a) (47) 


Thus the constants A and B are 



Hence the solution can be expressed in the following nondimens ional form: 


n^t rigt 

(w + wi)/(w elas t ic + = A'e + B*e 


(49) 


in which 


A' = jj(k - l)/rj[a/(l - aj - na 

/( n 1 - 02 ) 

b* - [( k - i)/ t 2 &A 1 - a 3 - n l 

1 (“1 - “2) 


(50) 


The ratio of the strain e at any distance z from the xy-plane 

(fig. 5) and. the elastic strain at z = 0, e n , can now be deter- 

u elastic 

mined from equation (12) together with equations (8a), (38), (44), and (49). 
Thus 


e / € Oelastic " - z ( aA / I )( 1 A.) ] [l/(l - a] (a 1 


njt 

e + 


B'e 2 J - 1 1 sin(jtx/L) + ( E-^E^fl - e ^ 2 ) + 


(V T l) + 1 


(51) 


in " Mch e °elastic " - p / AE l* 
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Equation (43) or (49) shows that for the column considered the deflec- 
tions due to linearly viscous oreep increase with time and approach infinity 
only as t approaches infinity for a = P/Pg < 1. Equation (49) is plot- 
ted in figure 7 for various values of P/Pg for a material whose elas- 
ticity and viscosity coefficients are given in the figure. The curves 
show that for P/Pg = 0.8 the deflections increase much more rapidly 

with time than for the smaller values of P/Pg . For the same material 

described in figure 7 , equation ( 51 ) is plotted In figure 8 for x = L/ 2, 
P/Fg = 0.5> and- z(aA/l) = ±0.1. The curves show that on the compression 

(concave) side of the beam the compressive strain increases monotonically 
with time, -while on the tension (convex) side the strain reaches a maxi- 
mum value in compression, whereupon it decreases monotonically through 
zero, finally becoming tensile. 

It may be noted that equation (49) reduces to the corresponding 
solution for a column constructed of a Maxwell element if Tg is taken 

infinitely large (see refs. 1 and 2). 


DISCUSSION 


The differential equation of bending of a beam or column was derived 
for a viscoelastic material whose mechanical behavior is analogous to 
that of the model shown in figure 2. This model, whose properties are 
defined by two elasticity and two viscosity coefficients, is the simplest 
model exhibiting instantaneous and retarded elasticity as well as pure 
flow. For materials whose properties are defined by more general linearly 
viscoelastic models (see fig. 1(a)) corresponding differential equations 
of bending can be derived in a manner similar to the derivation presented 
earlier in the text. 

The results of the analysis of the creep deflections of a visco- 
elastic beam under pure bending are shown in figure 3* It is seen from 
equations (8a) and (36) that the ratio of the deflections of the beam and 
the deflections of the corresponding purely elastic structure at any 
time is identical to the ratio of the total strain and the purely elastic 
strain of a bar under uniaxial tensile or compressive stresses. Consid- 
eration of the linear nature of the stress -strain-time relation (eq. (8)) 
also leads to the conclusion that for any bent bar the stresses are line- 
arly distributed across the depth since the strains have been so assumed. 
Hence for all cases of bending of a bar by static lateral loads the 
initial linear distribution of stresses across the depth of the beam 
remains -unchanged and can be determined by consideration of equilibrium 
only, even though the deflections increase continuously with time . 
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The analysis of the creep deflections of a simply supported column 
with an initial sinusoidal deviation from straightness resulted in the 
observations that the deflections tend to be infinitely large with time 
and that for the applied load less than the Euler load the deflections 
remain finite for finite time (see eq. (k9) and fig- 7) ■ Because of the 
linear nature of the general differential equation of bending (eq.. (20)) 
these observations are applicable to column problems similar to the one 
considered herein but with initial shapes other than the single sine 
wave. In the case of columns, as with beams, the stresses are linearly 
distributed across the depth (see eqs . (8) and (12)) . However, the 
actual distribution changes continuously with time since the applied 
moment which is proportional to the total deflection changes with time . 

The variation of strain with time at the midspan of the column at 
two points symmetrically located with respect to a plane normal to the 
principal plane of bending and containing the beam axis is shown in fig- 
ure 8. The compressive strain on the compression side of the beam 
increases monotonically with time; whereas the strain on the tension 
side reaches a maximum value in compression, then decreases monotonically 
through zero, finally becoming tensile. 

The deflection-time curves shown in figure 7 indicate that the 
deflections may increase very rapidly with time for end loads near the 
Euler load. It is interesting to note that the curves exhibit the pri- 
mary and secondary creep characteristics observed in experiments with 
aluminum columns reported in references 6 and 8. 


Polytechnic Institute of Brooklyn, 

Brooklyn, N. Y., April 18, 1952. 
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Figure 3 *- Creep strain for uniform tension or compression and creep 
deflections for pure bending. Ep = 5 X 10^ psi; E 2 = 13 X 10^ psi; 
Ap = 10^ lb-hr/sq. in.; Ag = 13 X 10 ? lb-hr/sq. in.; Tp = 2,000 hr; 

T2 = 100 hr. 
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Figure 8.- Creep Btrain at midspan of a simply supported initially curved 
column. P/P-g = 0.5* (See fig. 7 f° r material properties.) 


NACA- Langley - 1-14-54 - 1000 


